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SUMMARY 



In this paper the dpwnwash velocity distribution he- 
hind the wing was determined for the free jet and for the 
closed tunnel of "both circular and elliptic cross sections. 
The wing was placed at the center of the. tunnel. The theo- 
ry makes it possible to determine the downwash at any point 
in the jet. The computations were performed for points in 
the plane determined by the jet axis and the center-of- 
pressure line of the wing. The elliptic tunnel section 
chosen had an axis ratio equal to The downwash 

proved to be proportional to the wing lift and inversely 
proportional to the c ro ss- sect ional area of the tunnel,. 
Moreover, for the circular jet the downwash depends only 
on "the distance from the wing (see formula (29), figs. 8 
and 9) and for the elliptical jet it may be approximately 
represented by the product of a function depending only on 
the ratio wing span/jet width (fig. 25 for the free jet 
and fig. 35 for the closed tunnel), by a function depend- 
ing only on the distance from the wing (fig. 24 for the 
free jet and fig. 34 for the closed tunnel). 

The downwash velocities induced at the wing and due 
to the boundary effect of the jet have been computed for 
jets of different cross sections. The downwash velocities 
in the region b.ehind the wing have, on the contrary, re- 
ceived very little attention. A knowledge of these veloc- 
ities in the neighborhood of the wing, especially at the 
position of the tail surfaces is necessary, however, for 
the computation of the corrected pitching moments. At the 
same time it may be determined whether the curvature of 
the flow lines is sufficiently large to explain the exper- 
imental fact that the correction factor for the angle of 
attack is larger than that for the drag. 

. - - 

♦"Korrektur des Abwindes in Windkanalen mit kreisrunden 

Oder elliptischen Querschnitt en. " Luf t f ahrt f orschung , 
vol. 12, no. 8, December 25, 1935, pp. 250-264. 
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A, INTRODUCTION 



Up to the present the tunnel interference effect has 
"been determined only at the wing - more exactly, at the 
center of pressure of the middle wing section. The down- 
wash velocity at this point could he simply determined 
since it is exactly half of the velocity with infinite 
flow, and the latter may he presented as a two-dimensional 
potential flow problem: a system of images of the trailing 
vortices he ing set up at the boundary (references 1 to 4) . 
Figures 1 and 2 show these vortex images for open and 
closed circular jets* In the open tunnel the vortex images 
have the same sign as the original vortex; in the closed 
tunnel they have the opposite sign. 

The values of the downwash velocities from the wing 
outward to infinity wo ve no t known. For a rectangular tun- 
nel Glauert , who was interested in deriving a correction 
factor for the horizontal tail surface of an airplane, ob- 
tained a first approximation. For tunnels of other cross 
sections there was no comprehensive treatment.*) Seiferth 
has indeed tried to obtain better results for the circular 
free jet, but failed to obtain the proper correction fac- 
tors even at small distances from the wing. In the pres- 
ent paper, in part B, will be given the exact computation * 
of the induced velocities for circular jets, and in part C 
for jets of elliptical cross section. 



B. JETS OF CIRCULAR CROSS SECTION 
I. Statement of the Problem and Notation 



To a first approximation an airfoil may be replaced 
by a "horseshoe vortex," which corresponds to a uniform 
lift distribution (fig. 3).. Wings with nonuniform distri- 
bution may be replaced by several horseshoe vortices of 



Shortly after part B of our paper was completed, there 
appeared the Aerodynamic Theory of W. F. Durand , vol. II, 
in which Burgers considers the same problem but for a cir- 
cular jet only. The methods are entirely different. In 
the case of the free jet there is excellent agreement in 
the results. For the closed tunnel , however, ther'e are 
certain deviations which will later be considered. 
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different widths 2b (fig. 4). We shall therefore first 
carry out the computations for a horseshoe vortex. We 
shall use the following notation: The x axis coincides 
with the axis of the jet, the y axis lies along the span, 
and the z axis is perpendicular to "both and directed down- 
ward. Figure 5 shows the coordinate system used. The 
transverse part of the horseshoe extends from y = - b to 
y = + b , the trailing vortices run parallel to the x 
axis from x = 0 to x = ~ ». The circulation of the 
vortex will he denoted hy I\ The radius of the tunnel is 
r 0 . We set: 



(1) 



For infinitely extended flow the horseshoe vortex induces 
at the convex surface of the cylinder the tangential ve- 
locities u„ in the direction of the x axis and v+ 

perpendicular to u Q , and the normal velocity v nQ . 

In the case of the free jet the pressure must he con- 
stant at the jet boundaries; that is, the velocities u 0 
and vt 0 must vanish. We must therefore superimpose a 
velocity field which is regular within the cylinder and on 
the cylinder surface has the values -u 0 and " v t 0 » 

In the case of the closed tunnel the normal velocity 
v n 0 roust vanish. The superimposed field must therefore 

have the velocity -v no at the "boundary. 

II. Determination of the Downwash for the Free Jet 



We shall first consider the free jet. The boundary 
conditions for infinitely large -x may he satisfied by 
introducing the images of the horseshoe vortices outside 
the jet as shown in figure 6. This method is known from 
previous investigations (Prandtl, Tragf lugel theorie II, p. 
54), We denote the potential of the horseshoe vortex hy 
$ 0 , that of the image vortices by $ x . The potential 

($o + $i) will then satisfy the boundary conditions for 



y? + z 3 = r 3 = r 0 3 p 3 



~ = p sin <p;~ 



p cos cp 



= i and 



b _ 



= K 
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very large negative values of % = x/ r 0 # For f i hi ^ e val- 
ues of i there still remain tangential velocities at the 
cylinder surface. Figure 7 shows the velocities u 0 + u x 
and v't 0 vt . plotted against <p . The function u 0 + 

u x is even in <p and %\ the function v t(j + v ti is 
uneven in Cp and \. Both decrease rapidly with increas- 
ing ||| . To compensate these velocities, we introduce 
a potential $ 2 whose derivatives d$ 2 /dx and B$ 2 /r9cp 
vanish at infinity and on the convex cylinder surface; 
that is, for r = r 0 and p = 1, assumes the values: 

9$ 2 _ 3 (*o _±li 1 . H*2 = ~ 5 (Bp + J?i ) r 
9x " ~ dx * r9cp rdcp 

$ 2 must satisfy the potential equation A$ 3 = 0 which, 
in cylindrical coordinates, reads: 

i 9* 2 + 9 8 *a. + l 9!Sa. + = o (3) 

We shall first assume the function <£ 2 , instead of van- 
ishing as I = ~— to oe periodic in |. We shall 

then in our computations allow the period 21 = 21* r 0 to 
increase indefinitely. The assumed series for 

$ 2 = ^ t cos m cp | D mk sin J m (i ^ p) ( 4 ) 

where D m k are as yet undetermined coefficients and J m 
are Bessel functions of the first kind, satisfies the dif- 
ferential equation, since each term of the sum satisfies 
it. , 

In order that $ 2 may fulfill the "boundary condition, 
it is necessary that 



= /9^o + (5) 
V 9 x A-n 
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If we now develop ^ J ^ see fig * 7 ^ into 



Fo>urier series in Cp 



+ ( MoJLlA^ = _ r __ T I. CO cos . * (8) . 



o 

the boundary condition will read: 



-IE. 
4rr r 



.- Z cos m <p 2 D mk to cos ffj ^ Jm (i 
o m k l • \ ir y \ I / 

E l m (O cos m <p (5a) 



4tt r 0 m 

and "by comparing coefficients, we obtain: 

I Dmk f* -os I) J m (i = - l m (I) (6a) 

If we further expand the function l m ( £) into a cosine 
series in |, f — — is an even function in iy we 

again obtain by comparison the values of the unknown coef- 
ficients 



1 rr T l * ^ cos (H a ^ da 



* 



D 



mk 



I kn T (z krf\ , 

i* J ® V 1 - l 



In this way we finally obtain, by substituting (7) in (4) 
the expression for $ 2 



r B (3&i-p>iii(Sf. I s ) +l* . ■ 

f* (^) * -i. ^ ' /o 



Jr 

$ 2 = - -£- 2cos mcp E 

(8) 

If we now allow the period I to increase so that in 

the limit ^ = q; ^ = dq, then the expression for the po- 

t ent ial becomes : 

p ^sin q£ J m (iqp) , 1 + r °° n , x a /nN 

= "fir H cos mc P / q^J'Tiq - )""" dq TT 1 l m(cc)cos qada (9) 

mo 
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In the plane of the wing (that is, to a great. approxima- 
tion, in the xy plane for which 9 = 90° or 270°), 
there is obtained for the downwash velocity due to $ 2 • 

V r. . T sin q| J m (iqp) „ 1 t°° / \ 

"* do) 

and in the symmetry line of the wing, that is, for r = 0 



5$ !° sin oj | dq a +co 

r^ = * 4fr 0 -{ -~ jTTiO" n .£ ^ C ° S q a ** 

which is a relatively simple expression. 

The potential $ 2 was so determined that the condi- . 

tion: 

_ /B($ 0 + $j\ 
V3 X /r 0 " " v ax 

_ _ __T_ cos Cp r j i 

4nr 0 £ 2 +cos a <p L S W .P+1+K9-2K sin^ 



r o 



1 



} 



+ K 



{ 



J | 2 +1+K 2 +2K S in<p 

__— i — — + 

V | s + l+tc 2 -2K sincp 

1 } 



7 P +1+K 2 +2K 3in cp 
-2-smcp | ^-p--^ £ ~~ -~ -~ 

- -- ,) 

,/ £ 2 K 2 +1+K 3 +2K sincp v 

+ | 2 K 2 +l+K 2 "_2K sincp + 

y p k 2 +i+ k s +2 k s 1] (12) 
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was satisfied.. We must still show that the condition 

.,(«!il_!iL) (13) 
v d(f>; To \ r a cp / Tq 

is simultaneously satisfied. 

On the surface of the cylinder the tangential veloc- 
ity is 

r ? sin at 1 

( = t — £ m sin m<p / dq - / i (a) cos q a da 

Vrdcp/ 4TTr 0 m o 4 n J» m 

(14) 

This expression then must be set equal to the negative 

,3(« 0 + *A 

of the derivative (' — — ). For the derivative, 

^ J , we assumed the expansion 

/d($ 0 + $ x )\ r 

^ j r = Slmd' 1 ) cos m $ ( se e equation (6)) 

(14a) 

for simplicity 

l m (£»l) was written l m (£) 

Proa this it follows that the potential ($ 0 + $ x ) must 
have the form 

$ o + $ i = £7 2 T m(|.p) cos mcp (15) 
4TT m 



wljere 



lm(e.P) = - 9 ~i^ (16) 
,8(* 0 + $ x )> 



For the tangential velocity ( ) we then obtain 

\ r 3 cp / r 

by differentiation of (15) the equation 
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Now condition (13) requires that the derivatives (14) and 
(17) be equal and opposite. Since both are sine series 

in Cp, the coefficients of each harmonic must be equal to 

each other except for sign; that is, it is necessary that 

/ iijLAi dq t f l m (a.i) cos q a da = + r m U,i) (18) 

If we integrate the left side partially with respect to a, 
we obtain, using (16) the equation 

00 +00 

+ / sin q Uq ; / T m(a,l) sin q a da = T m (S,l) (19) 

and it is immediately seen that the left side is the Four- 
ier representation of the function T m (£,l); the poten- 
tial $ 2 therefore satisfies all condi t ions , 2 ) 

We shall briefly summarize what we have done so far: 
An airfoil in a free jet is compared with an airfoil in. an 
infinite flow, the potential of whose disturbing field is 
$i + $2* Due to tlie J et boundaries the additional down- 

a($ + $ ) 

ward velocities w%. = —i— l. £ - are induced, For the 

plane of the wing (that is, to a close approximation, in 
the xy plane) : 

a(», + 9 ) _ 3(* + » a ) ' 

— V = + r cf^ ((P = 90 t0 270 } (l9a > 

The derivative — £ is given by equation (10)" and — £ 

roqp r3$ 



It may be shown that the problem of determining a poten- 
tial having predetermined values for the tangential deriv- 
atives at the boundary reduces to the problem of determin- 
ing a potential whose values are given at the boundary sur- 
face; that is, to a boundary problem of the first kind, 
For this reason, both conditions (12) and (13) are equiva- 
lent to a single condition.. 
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Vr3<P/<p = 9 0 o Vdz / z =o 4 ^r 



1 

I 

V I 



Along the tunnel axis we finally obtain (see equation 
(ID): 



(20) 



+ * a ) = 

."a« 



Wv — 

• 4 ^ r o 



j e 2 * 2 + 1. 



L /l^IJj 4TTr 00 y J, (iq) TlJ. 1 

(21) 

For small values of £ the first term simplifies to 

\Bz / z=0 4nr 0 
y=o 

The downward induced velocities along the axis for 
different ratios of wing span tc jet diameter are deter- 
mined from equation (21). (See section IV.) 



III. Determination of the Induced Downward Velocities 

for the Closed Tunnel 



A similar procedure is applied to closed tunnels. Im- 
ages of the horseshoe vortices are first introduced with 
the sign of the circulation changed. (Pig. 2 shows a sec- 
tion through the tunnel at a great' distance from the wing.) 
Then the superimposed potential $ 2 * is determined in a 
similar way as was done for the open channel. We thus ob- 
tain 
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a* r o m r J m (iqp)sin q|dq it" 

— 00 

(22) 

in which J m « is the derivative of the Bessel function 
of the first kind and mth order. 



The functions g m are so determined that 

/a ttaS - = / B(^o + *iK 
A 8r A " V 3r I 



' r o 



coscp 



A — ( 



Ik 



1+K 3 _2K . sin9 V ! 2 K 3 +1+K 3 -2K sinCp 



e 



,/| 3 +l+K 2 -2K s incp' 



i ( 



1+K 2 +2K s in<p V | s K 2 +i+K 3 +2K S in<p 



£ 2 +cos 2 



y^ 3 +1+ K 2 + 2K sincp 7 . 

- 2 + sin cp (- T-- 

<P L V Tk 2 +1+k 3 -2K- S in(p 



| C0 3_|_ 

Lc» T Y V / 

V % K 2 +1+K 3 -2K 

L=± \ 

3 +l+K 2 +2K sine/ 



\/| 2k2 +1+K 2 -2K sinCp J t 2 K 2 +l+K S +2K sin?' 



+ sin cp ■ ~L ^ 

V I 2 +1+K 2 +2K ; sin9 y| 2 +l+ K2 -2K S in9 X 

+ K (-.- ----- - 1 - + - y -1 \]1 

V i +1+K 3 +2K sincp y| 2 +H-K S -2K sincp'JJ 

~-~ Z. g m U,l) cos racp (23) 



4TTr 



o fit 
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The induced velocity, arising from the tunnel "boundary- 
effect, in the plane z = 0, is 



4TTr~ 



P + 



/ Z 7 TT2 + / Z 7^ 



4rrr 



-foo 



y . m r J m (iqp)sin q £ aq 1 , . 
2, m sin mcpy ____ 7 __ — _ j g (a, l) sin qada 

(24) 



q J m ' (iq) 

and for the line of symmetry y = 0: 
r ]2 



+ 2K 



} 



- I*-' M-TrH" " Aja.D-m 1 a (25) 
4rrr 0o Ji'dq) tt ^ m 

From equation (25) are computed the induced velocities for 
various ratios of wing span to tunnel diameter. (See sec- 
tion IV.) 

IV. Numerical Results 



Sections II and III supply the theoretical basis' for 
the computation of the induced velocities due to boundary 
effect. These velocities and w^* are computed from 

equations (21) and (25) for any point on the axis of the 
wind tunnel. 

When the "blower velocity is V the change in the flow 
direction at each point is w k /V. The induced velocity at 
the wing will he denoted "by • The correction for the 

angle of attack of the wing itself is then 



r. 



(Xv = — ~ = 2 K 

k V 4TTr 0 V 
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circulation, K = 31. = ^iBg_span 



3r, 



tunnel diameter 



) 



or, substituting the lift coefficient 



ct k = 



w k, 



1 F 



8 F, 



(26) 



(F - area of wing, F 0 = cross-sectional area of jet) 



At a distance I from the wing the effect of the 
tunnel boundary on the downwash is larger than at the wing 
itself. The tail surfaces in that position would there- 
fore have a different effect on the pitching moment than 
would be the case in free flight. In the free jet the ef- 
fective angle of attack is reduced by the amount 



- Wt 



Pk = 



(27) 



in the closed tunnel it is increased by the amount 

Wjj.* - w k * 
Pk = f ~" (28) 

The upwash at the wing in the closed tunnel is numerically 
equal to the downwash of the same wing in the free jet, 
w k 0 * = ~ * V 

Figure 8 shows the ratio T = Pk/a^ plotted against 
l/2r 0 for the free jet. Figure 9 shows the similar rela- 
tions for the closed jet. } The computations showed that 
the quotients T and T* are practically independent of 
the ratio span/tunnel diameter. Since any distribution 
whatever could be built up by superposition of uniform 



Contrary to what was found by Burgers (reference 5) T*.for 
the closed tunnel approaches the asymptote from_above. 
W. L. Cowley and 0. A. McMillan (reference 6) have recent- 
ly improved the Glauert computation for the rectangular 
tunnel and have also shown that T* approaches the asymp- 
tote from above (for ratios of span to tunnel diameter less 
than 0,55). The x 4 curve the^e given would be similar for 
the case of a square tunnel but have its maximum at x/b = 
1 (b = width of tunnel). 
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distributions along spans of different lengths, the re- 
sults are also independent of the lift distribution. For 
the free jet, higher values are obtained than those ob- 
tained "by the Seiferth approximation. 



The induced velocity itself instead 
angle of attack, may easily he determined 
(figs, . 8 and .9) : 



Wk = i lo ° a U +T } 7 



* 

w k 



= - 1 

8 



JL 



(1 + 



T*) 



of the change in 
from the curves 



(29) 



It should he noted that by the above theory it is 
possible to compute the downwash at every point in the 
wind tunnel. Thus, for example, for a ratio of wing span/ 
free- jet diameter =0.45, the downwash was computed all 
along a parallel to the tunnel axis in the xy plane (fig. 
5) passing through the wing tip. The values differ by an 
insignificant amount from those along the tunnel axis. 

The downwash resulting from the jet boundary produces 
a curvature of the flow. We shall now investigate whether 
this curvature is sufficient to explain the experimental 
fact that the correction factor for the angle of attack is 
larger than the drag correction factor. 

Since, along the wing chord, the curves T and T 
(see equation (29) and figs. 8 and 9) may be replaced by 
their tangents f f (0) and f*'(0) through the origin, 
the downwash may be given by the expression 



w k = w ko _[i + f .(o) 2 -M 



(30) 



The curvature of the streamlines at the wing is therefore 

1 f'(0) v 'k Q 



R 



2r, 



(31) 



where f'(0) = 1.58 for the free jet and f*'(0) = 2.1 
for the closed tunnel. 



It is customary to write the correction for the angle 
of attack due to the boundary effect in the form 



a k = 



8 F, 



-a 
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For circular tunnels 6 = * !• This factor is computed on 
the assumption that the wing is replaced by a lift line. 
Since the wing, however, has finite chord, the jet curva- 
ture does affect its. camber characteristics. In a free 
jet the camber is decreased - that is, the effective angle 
of attack is decreased; in a closed tunnel it is increased, 
and therefore the effective angle of attack is increased. 
Seiferth (reference l) investigated the effect of jet cur- 
vature on a profile having the form of an arc of a circle. 
If t denotes the profile chord, then the change in the 
angle of attack of a wing section due to the streamline 
curvature is of the amount 

Aa=-| (32) 

If we substitute the value for E given "by equation (31) 
we obtain for the angle of attack the correction 

t £t(0) w k Q t f'(0) , v 

. Att = 2r~ "4- "V- " 2r7 ~4~ ( 33 > 

The total correction for the angle of attack for a free 
jet is then 



+ A a = 



1 + 



that is, the factor |s| increases by 

100 (£) G^)T'-»' (35, 

If we compare this result with those of Knight and 
Harris (reference 7), there is found to be satisfactory 
agreement for the ratios wing span/tunnel diameter * 0.45 
and 0.6; for larger ratios of K = 2b/2r Q , the values ob- 
tained are too small. We must consider, however, that our 
computation was made only for the tunnel axis, so that no 
account was taken of the variation along the span. For 
large spans, for example, K = 2b/2r 0 - 0.3, this varia- 
tion along the span would increase the factor 6 by 9 
percent even if no account were taken of the curvature of 
the jet. Moreover, the test results were obtained by cx- 
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trapolation to infinitely large ratio of jet area to wing 
area and are therefore not very accurate. 



C. JETS 07 ELLIPTIC CROSS SECTION 
I. Statement of the Problem and Notation 



Due to the fact that they make the testing of large 
models more convenient, free jets of elliptic cross sec- 
tion are now preferred to circular jets. For uniform lift 
distribution, the induced velocities at the wing arising 
from the boundary effect were computed by Sanuki and Tani; 
Rosenhead made the computation for elliptical lift distri- 
bution (references 3 and 4). We shall determine the in- 
duced velocities behind the wing in the same way as was 
done for the circular jet« 

We again first replace the airfoil by a horseshoe vor- 
tex, since nonuniform lift distributions may be built up 
by a combination of horseshoe vortices. 

The direction of the coordinate axes and the notation 
for the various magnitudes of the horseshoe vortex are the 
same as in section B. The elliptic jet has its large a-^is 
equal to 2a and its small axis to 2bi. (See fig. 10.) 
The eccentricity is denoted by c (c = J a 2 - b L 2 ) . With 

infinite flow there are induced at the convex surface of 
the elliptical cylinder the tangential velocities u 0 in 
the direction of the x axis and vt 0 perpendicular to u 0 ., 

and the normal velocity v n 0 • 

In the case of the free jet the pressure must be con- 
stant at the jet boundary; that is, the velocities u 0 and 

v t must be made to vanish. We must therefore superimpose 
0 

a velocity field, that is regular within the elliptic cyl- 
inder and on the boundary has the. values ~u 0 and ~ v t 0 * 

3?or the closed tunnel it is necessary that the normal 
velocity v nQ vanish* The superimposed velocity field 

must therefore have the value ~v nQ on the surface of the 
elliptic cylinder. 
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II, Downwash in the Free Jet 



1 . T h e o r e t i c a 1 cons i d ,e r a t i o ri_ s_. - We shall proceed in a 
similar manner as for the circular free jet*. In that case 
we first set up the images of the horseshoe vortex at the 
cylindrical ."boundary and were thus enabled to satisfy the 
boundary condition for infinitely large negative values 
of x. In the case of the elliptical jet, we must similar- 
ly arrange a system of vortex filaments so that the bound- 
ary condition at infinity may be satisfied. Since the ve- 
locity u 0 approaches zero at infinity whereas the other 
component of the tangential velocity v^ 0 remains finite, 

it is necessary to determine the vortex filaments outside 
the cylinder ir t such a manner that they make vanish 

at infinity. ' 

If we consider a section of the jet at an infinite 
distance behind the wing (fig. 11), we have an ellipse 
with two single vortices. We must now set up a system of 
vortices in the space outside the ellipse in such a manner 
that the velocity induced by these vortices and by the vor- 
tex pair within the ellipse is always normal to the bound- 
ary. To determine the position of these vortices we trans- 
form the ellipse into a rectangle. . In figure 12 B and B 1 , 
denote the foci of the ellipse. Let there be a vortex at 
P x . By means of the function 

y + iz = c cosh {r\ •+ it) (36) 

the ellipse is transformed into a rectangle (hatched area 
in fig. 13). The segment BC = fiu goes over into segment 
B*c* « t| 0 on the r\ axis. The upper half of the ellipse 
is transformed into the rectangle B*C*C , *B t * lying above 
the r| axis; the lower half into the rectangle B*C*C ! *B f * 
lying below the r\ axis. The side of the rectangle C'*D* 
G*D l *cT I * corresponds to the boundary of the ellipse. In 
order that only normal velocities arise at the boundary 
C T *D* C*D'* C 1 * and that the flow in the yz plane pass 
through C'B^B without disturbance, it is necessary to 
form images as indicated in figure 13. Let the complex 
coordinate of the point P x * be p, x . The image vortices 
then lie at 
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1) p, x + i2nTT + 412^0^ 

2) (2t) 0 - + iSnTf + 4mq 0 

3) - M-i + i2nrr + 4rnn 0 

4) (2tiq + ili) + i2nrr + 4mr) 0 



(37) 



where jj. x denotes the conjugate of \x x . Assuming that 
Pi* = jJix denotes" the position of one trailing part of the 

horseshoe vortex, the other vortex will lie at point p 2 * 
whose coordinate is |Jl 2 = in + jl x . The position of the 

corresponding image vortex may he' found with the aid of 
equations (37) by substituting in + [i x for \x x and 
-in +, |jLi f or pL x . 

In what follows we shall consider only the case where 
the airfoil lies in the xy plane; then the z coordinate of 
the trailing vortex is zero. If the span 2b is now great- 
er or equal to the distance between the foci, then in the 
r\l plane all vortices will lie on the straight lines £ = 
nn; if the span 2b is smaller than the distance between 
foci, the vortices will lie on the straight lines r\ = 2m 
r| 0 , where m and n run through the whole range of in- 
tegers from -co to +oo. We must still determine the po- 
sition of these vortices in the yz plane. Since y = c 
cosh n cos £ and z = c sinh r| sin £ (see equation (36)), 
all vortices will lie on the y axis for the range 2b 
= 2c, as £ = nn; and for 2b < 2c, all vortices will 
lie on lines for which £ = constant 4 nn; that is, on 
confocal hyperbolas. As 2b approaches zero, all vortices 
lie on the z axis. Thus we have determined the positions 
the vortex filaments outside the boundary must have in or- 
der to satisfy the boundary condition for x — s> - co. 

Now let $ 0 be the potential of the horseshoe vortex, 
$ x the potential of the vortex filaments in the outer 
space. Figure 14 shows the vortex images for the case 
where the span is equal to or greater than the distance be- 
tween foci. The figures indicated in the diagram arc to 
be used later for working out an example. In the zy plane 
it is necessary to prolong the vortex filaments in some 
manner since a vortex filament must have no free ends and 
since, moreover, they cannot be prolonged parallel to the 
x axis because in that case the boundary condition for 
x +oo would be impossible of fulfillment. For the case 
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where the span is smaller than the distance between foci, 
the position of the vortex filaments will be as shown in 
figure 15. In this case.it is simpler to unite the fila- 
ments rather than prolong them. This combining of the 
filaments would have, been impossible in the first case 
since the transverse part would have cut the ellipse and 
so would not have lain in the outside region. In the case 
of infinitely small spans vortex pairs are obtained in 
the space outside the elliptic jet (fig. 16) ♦ 

The potential ($ 0 .+ $i) satisfies the boundary con- 
dition at infinity. In order to reduce to zero the ve- 
locities still remaining (u 0 + u x ) and (vt 0 +■ v t x ) » it 

is necessary to superimpose a potential $ 2 . Since $ 2 

must satisfy the above conditions on the surface of the 
elliptic cylinder, we introduce elliptical coordinates 
for greater convenience in treating the problem: 



x = e 

y - c cosh r| cos £ 
z = c sinh n sin t 



(38) 



In any plane .. x = constant the curves r\ = constant rep- 
resent confocal ellipses and the curves £ = constant, 
confocal hyperbolas. On thfe cylinder surface r\ = con- 
stant = n 0 • The potential equation which $ 2 must satis- 
fy then reads, in elliptical coordinates: 

£V' ~ ^!M.'*tif\ m0 (39) 

d£ 2 c 2 (co S h 3 n-cos 2 o v a-n 3 U '- 

In addition, the "boundary conditions require that 

, . • ^0 ^0 

and 

(a"A 0 .- - (-- 3 -r-), 0 <«> 

Instead of allowing $ 2 to vanish i — oo , we 
shall first set up $ 2 as a periodic function in \ and 
then in our computation allow the period 21 = 2 1* a to 
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grow indefinitely. To solve the differential equation 
(39), we assume the product: 

$ 3 = - {- Z 2 C Uik sin (If l)M Utk (» M Vfk <T)) (42). 

The function M depends only on £, the function M on- 
ly on n* Substituting one member of the above equation 
in differential equation (39) , 

- sin 1") M„ k U v k +— -3- 1 - — sin (?3E A 

I 2 \l J v ' ]c c 2 (cosh 2 n-cos 2 C) ^ l ^ 

S, k (Ti) + ^1— - M u>k (£)) = 0 (42a) 

it is seen which equation the functions M and M must 
satisfy: 

— — |-ar M Vtk (t)) + M Uik (£) 

2 2 

= k__TT_ c s (cosh 2 ^ C0S 3 £ ) Mu(k (J.) (t,) 

If 

= ^-g 11 - |- (cosh 2ry-cos 2 0 M pik (0 M Ufk {r\) (42b) 
I 

Changing the order of the terms and dividing through by 
M M, we obtain: 

^ + — " 2 — -o- cos 2 J 



d ? 2 M v>k ({) v 2 2 

d2 M„ k (ti) 1 k 3 tt 3 c 3 



d Tl 2 M U(k (t)) i 3 2 



(42c) 



Since one side depends only on the variable £ and the 
other side only on r\ , the equality of the two sides means 
that they are constant. We denote the constant by -V, 
and thus obtain the equations: 
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3 — + (~--- ~ cos 2$ + v) M u>k ({) = 0 (43) 



d2 M V,k (Tj) A 2 TT 2 C 2 



dTl 2 ^ l : 



— cosh 2^ + M P)k (r))= 0 (44) 



If = iT is substituted in the second equation, it as- . 
simes the same form as the first* The solutions of the 
first equation for real values of £ are therefore, «t>n 
substitution of a pure imaginary variable, solutions of 
the second. These differential equations were first in- 
vestigated by Mathieu, and their solutions are known as 
Mathieu functions (references 8 and 9). For our problem 
we need consider only such Mathieu functions as. are periodic 
in £. Periodic solutions, however, occur only when there ex- 
ist* definite relationships between the parameter v and 

2 2 2 

tLaf , s0 that i? is determined when - T1 >-- %- "is 
given. 

In addition to satisfying the differential equation, 
the potential $ 3 • must also satisfy the boundary condi- 
tions (40) and (41) . Now the problem of determining a po- 
tential that should have, given tangential derivatives at 
the boundary, is equivalent to the problem of determining 
a potential having pre-assigned values on the "boundary; 
that is, to a boundary problem of the first kind. : There- 
fore, both conditions (40) and (41) are equivalent to one 
single condition and it is therefore sufficient that the 
condition: 



(ae% 0 --( — iT V *"* ,Ce, > , « ) 



(45). 



is satisfied, and therefore also 



The constants C v ^ must be so determined that the bound- 
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ary condition (46) is satisfied. Since the Mathieu func- 
tions i^ f k ( 0 form an orthogonal system, similar to the 

harmonic functions, to which they reduce when the eccen- 
tricity of the ellipse approaches zero, the following re- 
lations hold true: 

2TT 

/ V>k ( £> M v m ,k (0 *C = 0 for n 4 m 



2TT 



/ M % n ,k CO dC = Tt for n > 0 

o 

2TT 

and / M 2 Vo , k ( £) d£ = 2tt 

We shall therefore multiply equation (46) "by a Mathieu 
function My t ( 0 and - integrate over the region 0 ^ 
C = 2rr. 

- f z c v k cos flyL £| ¥ k (r, 0 ) 

4tt k v ' £ l a \l a/ 0 

= - %Fa £ ' 7 ( ^o'0 CO d£ 



(47) 



(47a) 



whidi , divided "by t~ » gives: 

4rra 



? °v,k 7^ cos 



_ 1 



2TT 



*■ / » (e.ti 0 .c) M v>k (o d^ (47D) 

0 

We shall now make use of the orthogonal properties of the 
trigonometric functions "by multiplying the equation "by 

cos ( l ~L IN an< i f integrating over the region - \* = 1 = + 

^ \ I a/ a 
1 . We thus obtain: 



+1 



(•Ho) = 



n +1 * t 2TT * 

*■ 4 tcos (f£ |)£ / Ke.Ti 0 .£) M U(k (o da d (l) (47c) 
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or ' c i;,k = 

* 

(48) 

Substituting these coefficients into equation (42) and 

I * 
setting - = a and £ = £ under the integral signs to 

avoid confusion, there results for the potential $ 2 the 
expression: 



(f* I) 



+!> v -, 2tt 



£ / + [cos (^ a )da^ / Ka.Ti 0 . £*)M v>k ( £*) d £*] 



M l?,k (0 M u>k (ti) (49) 
rioc 

limit J?- = q, —jr = dq , and thus obtain for the poten- 
tial the final expression: 



We now allow the period \ to increase so that in the 



3 4TT u 

oo sin (a — 1 + 00 , 2tt 



" 7 ^TjT) / cos *ada-/ F(a,n 0 ,£*)M (£*)d£* 

The Kathieu functions (see (43) and (44)) will then satis- 
fy the differential equations: 

df 3 "~ + C q8 i^ C0S 2 * + V W» = 0 (51 > 

and 
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' U2 2P °° Sh 2X] + V) M V,q (T 1 ) = 0 (52 > 

The downwash velocity required is now given by 

8$ _ d$i 3^2 _ 3$^ cosh ^ sinj 9*a. 

6z 3z 3z 8z c(cosh r]-cos £) 3 r\ 



+ ___sinh_U_cos_l__ ||a (53) 
c ( cosh r)-cos £ ) a s 



For points in the xy plane where the y coordinate is small- 
er than the distance from the foci to the origin, that is, 

- c ^ y 5? + c so that r\ = 0 , 

5* = 3*1 + 1 ( 54) 

82 a z c sin i 3ti 

and for points whose y coordinate is greaterthan this dis- 
tance, that is, -c = 7 = •+ c so that £ = 0 

9$ 3*! 1 . , , 

— = + — ~ ( 55) 

3z 3z c sinh r\ o £ 



N^§.r i c a^_e xampl e . - An elliptical jet section was 

chosen with axis ratio^T a/b x = On the jet boundary 

the elliptic c«o;o rdinate r\ has the value n 0 given by 
the relation: 

~ = coth n 0 = (55a) 

The computations were carried out for the two following 
cases : 

1) Wing span equal to the distance between foci. 

2) Very small wing spans. 

3($ Q + $ ) 

We first require the tangential velocity ~r — — • 

For the first case the vortex filaments outside the jet 



4) 

The notation b x was used for the minor axis of the el- 
lipse to avoid confusion with the wing span. 
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all lie in the xy plane. The velocity at the jet "bound- 
ary (r\ •= T| 0 ) is therefore, 



4 sin I 



cos £ + -5 
a 



c U 

cose- - 
a 




where 



n=oo 
+2 2 
n= l 



-2 + — 



cos £ + d n 



(i) + Gr) sin2 ^ (cos C + a n ) : 



cos £ - d n 



0 + Sin2 ^ + (C0S * - d ^J, 



(55) 



co sh 2xi T) 0 
^ n " cosh T) 0 



n = 1) 2 , 3 • • 



(56a) 



For the second case we have pairs of vortices in the xz 
plane (see fig. 16); the velocity at the jet boundary is 
given "by 

aC^+fci) _-F2V n=oo |^ sin £ -f n J cosh 2n r\ 0 



4na 2 n= - 



=S (57) 
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where 



sinh 2n r\o ^ , 



The effect of the vortex filaments lying outside the 

T O f rl o^-^ao aO« *tt T +■ *U 4 V\ Vt n n { vt m J j » WM . -C* > 4- "U — _ X. T X, 

is sufficient for practical purposes to consider, out of 
the infinite number of vortex filaments lying outside the 
jet, only the first three'. 

The velocity on the | direction is an uneven func- 
tion of £ . We must choose the corresponding Mathieu 
functions M Vi .q(£) to satisfy the differential equation 

(51). S, Goldstein (9) has computed several Mathieu func- 
tions using a notation for the constants different from 
ours. To avoid confusion, we shall employ the index G 
with the Goldstein parameter: 

c 2 

The uneven periodic Mathieu functions Goldstein denotes by 
s £2n-fl (0 since they reduce to the harmonic functions 
sin (2n + l) £ when the eccentricity of the ellipse ap- 
proaches zero. 

The Mathieu functions may he developed into Fourier 
series : 

■ e 2n+l <0 = 2 B 2n+lf2m+1 (~) (2m + l) £ (59) 

The function s e x may he taken from the tables computed 

by Goldstein, noting that q & is negative for real values 

of q (reference 9, p. 304), The functions s e 3 , s e 5 , 

etc., as far as we know, have not yet been computed. Using 
Goldstein's method, we first computed. the values of v 
and then the coefficients of the Fourier series to obtain 
the functions s e 3 and s e. s . 5 ) Figures 17, 18, and 19 

show the variation of the Fourier coefficients of s e x , 
s e 3 , and s e 5 with qg. 

We require, in addition, the corresponding Mathieu 

functions with imaginary arguments. These will be denoted 

_ _ _ , 

The computations were carried out by P, Riegels, In com- 
puting the values of v the work of E. L. Inco (reference 
10) was used. 



26 



N.A.C.A. Technical Memorandum No, 801 



by the ietters Sgn-fl (r\) • B. Sieger (reference 11) has 

shown that between the corresponding functions there ex- 
ist relations such that 



s e QC (£) d£ (60) 

2n+l.V 



where 

r cos $ = c cosh n cois £ 
and r sin £ = c sinh n sin g 



(61) 



With the aid of an addition theorem for Bessel functions 

( reference 12) , the integration of (60) may he performed. 
The theorem is: 

X U r U (V) m=o\ ( _ l} m Ap^ 

C\ (- cos 2 0 (62) 

where P(u) denotes the gamma function and C° m (p.) the 

coefficients of cc m in the development of (l - 2a p + a 



a 3 ) U according to increasing powers of a. Prom equa- 
tion (60) there results: 



SSn+l.qOn) = sinh <n "z* h 2n +l ,m( l) *m (| | e ' n ) I m(|- J e"^) 



m= 1 

(63) 



with 

m."f* 1 

b 2n+l,m = m k B 2n+l,l~ B 2n+1 , 3 • V B 2n+l-.2m-l^ ( 64 > ■ 

For the function the series is convergent. Tor higher 

orders it is often necessary, for values of r\ that are 
not too small, to obtain the function directly by integra- 
tion of (60), since the accurate computation of the higher 



order Bessel functions I m ( - o e 7 ! ) would take too much 
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time. We shall limit ourselves in our computations, to 
the determination of the downwash velocity behind the 
wing. Tor this we have the value (see equation 54)): 

» « »tL + L_ T |i«L ( 65 ) 
dz 3z c sin £ B r\ 

Prom equation (50), substituting 

V (0 = S 62n+1 ^ ( " (66) 
. . M U,q ('O) = §-2n+l,q li)) 

there results: 



3 2 4TT ii TT Q 



iin (q I) i.' 2 n+l,q(0) • , 

— / j cos qada 

1 S2n+l.q (t W -» 1 

/ F(o,Ti 0 .C*) ■ e 2n+1 . >q (S*)d^se 2n+lfq (Oj dq (67) 
We therefore require to determine only the quotients: 

m 2n+l,q = s rirl" C ' 

This quotient is plotted in figures 20 and 21 for n = 0 
and n = 1 . 

In the computation it turned out that the first term 
of the series in n contributed the esse ntial part. 
This is due to the fact that for small values of q only 
the first and third terms are significant in the Mathieu 
function expansion s &2Ti+l,q (£)• 

3« Re s ults .- With a free elliptical jet with J 2 ax- 
is ratio, the downwash velocity due to the jet boundary 
was determined along the jet axis for the case of a wing 
of very short span and one whose span ttas equal to the 
distance between the foci. It was assumed that the wing 
had uniform lift distribution. For. the finite wing there 
was, in addition, determined the mean value of the veloc- 
ity behind the wing: 
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/ — z~ d y 

^ = -c *_ z _ 



(69) 



as a function of the distance l/2a from the center of 
pressure line of the wing (fig. 22)". 

It is usual to represent the mean downwash velocity 
in the form: 

= — 6 (70) 
V 8 F 0 

where c a is the lift coefficient of the wing, F Q the 

cross-sectional area of the tunnel, V the blower veloc- 
ity, 6 . a nondimens ional coefficient which is a function 
of the distance I from the center of pressure line of 
the wing and of the ratio wing span/jet width. Figure 
23 shows 5 as a function of the ratio wing span/ jet 
width. For comparison, the corresponding values for a 
circular jet are also indicated. If the coefficient 6 
in equation (70) is split up by using 5 0 to denote the 
value of 6 at the wing: 

^= °f ~8 0 (1 + T) (71) 

it turns out that the curves T (l/2a) practically coin- 
cide for b = 0 and b = c, so that it may he assumed 
that T is nearly independent of the span for one and the 
same jet section. It is therefore possible to determine 
the downwash velocities for^all spans between 0 and 2c 
from the foregoing computations, provided 5 0 is known, 

Figure 24 shows T. For comparison the corresponding val- 
ues for the circular jet of equal width (2a = 2r Q ) are 

also indicated on the figure; in that case, 8 0 = !• The 

value of 5 0 for the elliptical jet is known'from the 
computations od Sanuki and Tani (fig. 25). 

In order to distinguish the effect of the elliptical 
cross section on the downwash behind the wing we shall 
compare with a circular jet of equal cross-sectional area. 
For a wing of very small span we may obtain the value 6 
of the elliptical jet approximately by taking the value 
for a circular jet of equal area and considering that §(0) 
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is not equal to 1 but is 6 0 . Figure 26 gives the result 
and shows that the approximation is good, which result 
may he expected when the eccentricity of the ellipse is 
not too large. 

. . j x n 4- ~ * J? * V» ^ w-S v\ p»» •? 4- n o r\ -f i nt. P V— 

111 IUt?cli3 UI XU£ U U.O lUUiUCiil/9 V J- vj.**S» »« -»-**£> — " » w - 

est to know the distribution of the downwash along the 
span. Figure 27 shows the variation of the downwash along 
the wing span having uniform lift distribution. 

III. Induced Velocity in Closed Tunnel 



l2._The0.ret ical_con si derations.- In an infinite flow 
the airfoil would induce at the surface of the elliptic 
cylinder, which corresponds to the tunnel boundary, the 
normal velocity v nQ . Since the normal velocities at the 

cylinder surface must vanish at the boundary, as the air 
can flow only along the wall, we must superimpose a veloc- 
ity field which has a normal velocity component ~ v n 0 at 
the surface. We proceed in the same way as was done for 
the elliptical free jet. We first superimpose on the po- 
tential field of the wing $ 0 a field 9 1 which is pro- 
duced by the so-called "vortex images" and which assures 
that the boundary conditions are satisfied at infinity. 
The position of these vortex images is the same as for the 
free" jet but the sense of rotation is partly different.^ 
To determine the latter we consider a section through the 
tunnel at an infinite distance behind the wing (fig- 28). 
The vortices outside the ellipse must, together the 
vortex pair within the ellipse, induce on the ellipse/a 
velocity which is always in the direction of the tangent 
to the ellipse; that is, which has no component normal to 
the boundary. We may see these relations clearer if we 
transform the ellipse into a rectangle (fig. 29). The 
curved boundary of the ellipse then goes over into the 
straight line 0**0'* • The flow must now take place along 
this straight wall. 1% may be seen that, contrary to Vrhat 
was the case for the free jet, the sense of rotation must 
change along the horizontal, lines but remains the same 
along the vertical lines. Figure 30 shows the position of 
the vortex, images in the xyz system of coordinates for 
the case where the span is equal to the distance between 
the foci 2c. 5 ) The sense of rotation is indicated. For 



In this case the point 
origin* 



? x * in figure 29 moves to the 
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spans 2b < 2c and 2b = 0, the corresponding figures 
may easily he obtained by comparing with figures 15 and 16< 

The potential ($ 0 satisfies the boundary con- 

dition at infinity. To reduce the remaining velocities 
( v n 0 + v nJ to zer0 i* is still necessary to superimpose 

another potential $ 2 , which is determined by the condi- 
tion that the equation 

1¥=" d¥*80-- "n 0 + v ni) (73) 

holds true at the cylinder surface.. For the potential 
we obtain in a similar way as for the free jet: 



u 0 



r~= 7~r~ / J sin q a da 



n 0 



2tt 



£ / F* (a,n 0 ,n M u>q (£*) d£*}l 

o J 

where | , r) , £ denote elliptic coordinates. On the cyl- 
inder surface n = constant = r\ 0 . The coordinates \ and 

£ appear in the inner integrals as variables of integra- 
tion. To avoid confusion the integration variable - is re- 
placed by £* and g/a by a. M Vtq U) and M pjq (r]) 

denote the Mathieu functions with real and variable argu- 
ments which satisfy the differential equations (51) and 
(52). The function F* (a,T) 0 i£*) is connected with (v nQ 

+ v n x ) tlie relation: 

F* (a.Tlot-C ) = ^ Vcosh. no-cos 3 £ (74) 

4na 

If the components v y and v 2 of the sum are introduced 
in place of the sum, there results: 
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r*(cc,Tio. ~ { v y sinh 'Ho cos 5*+ v 2 cosh n 0 sin J* J 

4TTa (75) 

The induced velocity we are seeking is, in general, 
given by equation (53), for points in the xy plane by equa- 
tions (54) and (55). We shall limit ourselves to the case 
where the span of the wing does not exceed the distance 

between the foci. We then require — 8-' For this purpose 

the value of ^-^ for ri = 0 shall he specified and out 
of) * 

of the solutions of the Mathieu differential equation, those 
suited to this problem chosen and substituted: 

(a „ ) - ST 5 „- / S~~— "J Sl ° (* a) / q a da 

, 2TT 

J / I*(a.Ti 0 .C*) 3 e 2n+1 (£*) d£*] 

0 

s e 2n+l,q (0} (76) 



The quotient 



m 2n+l,q = gt 



a 2n+l f q (t) 0 ) 

for n = 0 is shown on figure 31 plotted for the case 
where the axis ratio is ,/lT. 

2. Results^ .- A jet section was chosen whose axis ra- 
tio was a/b x = J 2. The induced velocity along the tunnel 
axis was determined for two cases, namely: 

1) Wing span equal to distance between foci. 

2) Very small wing spans. 

The lift distribution was uniform. For the wing of finite 

span, there was, in addition, determined the mean velocity 

behind the wing as a function of the distance l/2a from 

the center of pressure line of the wing (fig. 32). The 

downwash velocity for the closed tunnel is negative, that 

_ 

'The computations were actually carried out by F. Riegels, 
for whose assistance I take this opportunity to express my 
thanks. -* 
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is, it is really an upwash velocity, which leads to an in- 
crease in the angle of attack of the wing. For a given 
lift in a closed wind tunnel the angle of attack measured 
is therefore too small* 

The mean downwash velocity in a closed tunnel may 
again he given in the form: 



c a 3 1 « 

T ? o 5 



(78) 



Figure 33 shows 8 as a function of the ratio wing span/ 
jet width. For comparison the value of 8 for a circular 
jet of equal width is also indicated. 



If in equation (78) the coefficient 
denoting its value at the wing by 8 0 > 



V 



a ' o 



(1 + T ) 



8 is split up, 



(79) 



the result is, as in the case of the free jet, that the 
curves T (l/2a) for b = 0 and b = 2c run sg close 
to each other that for all practical purposes T may "be 
assumed independent of the span for one and the same cross 
section. This means that in this case too it is possible 
to obtain the downwash velocity for all spans between 0 
and 2c provided only 8 0 is known. Figure 34 shows T*; 
for comparison, the corresponding value for a circular sec- 
tion of equal width is also shown. The correction factor 
8 Q for uniform lift distribution was computed by Sanuki 

and Tani (reference 3); it is shown in figure 35 as a 
function of the ratio wing span/tunnel width. 

At the distance I behind the wing the induced veloc- 
ity produced by the tunnel boundary is greater than at the 
wing itself • A horizontal tail surface at that position 
therefore has a different effect on the pitching moment, 
at a given angle of attack, than it would have in free 
flight. In the closed tunnel the effective angle of attack 
is increased in the ratio 



A w k * 



" 8 F 0 60 



(80) 



Translation by S. Reiss,. 
National Advisory Committee 
for Aeronautics. 
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Figs. 1,2,3,4,5,5,7 
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Figure 1*- Vortex Figure 2.- Vortex Figure 3.- Wing Figure 4.- An air- 
images images with foil 
for the open for the closed uniform lift dis- "built up of horse- 
circular jet. circular tunnel. tribution (horse- shoe vortices. 

shoe vortex) 



* » Figure 5.- System of 

coordinates 





Figure 6.- 
Images of 
horseshoe 
vortex 
at the 
„ boundary 

180° ~ . , 

of the 
free jet 
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figure 8.- Value of t along Figure 9.- Value o± r* along jet 
the free jet axis* axis in closed tunnel. 




Figure 10.- Coordinate system. 
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Figure 11.- Section through 

tunnel infinitely 
distant "behind the wing. Trail- 
ing vortices are indicated. 
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Figure 12. Figure 13. 

Figures 12 and 13.- Transformation of ellipse into rectangle. 



1I.A.C.A. Technical memorandum No. 301 Figs. 14,15,16,17,18,19,20,21 

Figure 14,- Images of horseshoe vortex at 
the elliptical jet bound- 
ary. Span greater than or equal 
to the distance "between foci. 





Figure 15.- Images of horse- 
shoe vortex at the 
elliptical jet "boundary. Span 
smaller than distance between 
foci. 




Figure 16.- Images of horse-^^^ 

shoe vortex at f 
the elliptical jet boundary^/ 



for infinitely small span 
of wing. 




Figure 18. 



4 - 5 —'Cddst?** 0 



Figures. 17, 18, 19. -Coefficients 

of the period- 
ic Mathieu functions se^,se^,se 5 
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Figures 20,21, 
Quotients of 
the Mathieu 
function. 
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Figs. 22,23,24,25,26,27,31,32 
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Figure 22 .- Downwash "behind airfoil in 
elliptic free jet. 
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Figure 23. -Variation of the factor 6 

with ratio wing span/ jet 
width and with shape ^f cross section. 
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Figure 24.- For the confutation 

of the mean value of 
the downwash w^/V "behind an air- Figure 25.- Correction factor 6 0 for 
foil in an elliptic jet(l/2a= the computation of thedown- 

distance from pressure line wing) wash at the wing in elliptic free jet. 
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Figure 31.- Quotient of the Mathieu 

function ^plotted against q. 
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Figure 26.- Effect of shape of jet. 
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Figure 27.- Distribution of down- Figure 32.-" Downwash" behind airfoil in 
wash along wing span, closed elliptical tunnel. 
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Figures 28,29.. 

Transfor- 
mation of 

ellipse into 
rectangle. 
Images for- 
mation in 
closed tunnel. 
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Figure 28 



Figure 29 




Figure 30.- Images of horseshoe vor- 
tex at thehoundrjr of the 
closed tunnel of elliptic cross sec- 
tion. Span equal to distance "between 
foci # 



Figure u5.- Correction factor 

6 0 for computation 
of "downwash" at the wing in 
closed elliptic tunnel. 
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Figure 34.- For computing the mean 
value of the "downwash" 
wjv "behind an airfoil in closed 
tunnel of elliptic cross section( l/2a~ 
distance from pressure line of wing)# 
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Figure 33.- Variation of factor 
6 with ratio wing 
span/width of jet and with shape 
of jet section* 



